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The Dzialoshinskii model of periodic and helicoidal structures has been analyzed without 
neglecting of the amplitude function oscillations. The amplitude function oscillations are shown to 
be important for understanding of the nature of the phase function. Analytic consideration is carried 
out in the limit of small anisotropy (neglecting the cosine term in the Hamiltonian). Surprisingly, 
the phase jumps survive even in the limit of the vanishing anisotropy. 

I. INTRODUCTION 

The theory of modulated structures in macroscopic ferromagnets [l|, ferroelectrics [3| and metal alloys [3| has a 
rich history. The number of works dedicated to a study of the modulated structures in small samples is significantly 
less [1]. The standard approach to the theory of the modulated structures (including the incommensurate ones) uses 
the constant order parameter amplitude function approximation Q, 0|- Our consideration of the bounded sample 
(thin film, in particular) requires us to drop this approximation: the order parameter must satisfy some boundary 
conditions and it is unlikely to find a solution with constant order parameter satisfying these conditions. That is why 
we have to consider a generic case of connected set of the nonlinear equations for the order parameter amplitude and 
phase functions. We will show here that the solutions of this system are significantly more complicated, than in the 
constant order parameter approximation. 

II. FREE ENERGY AND THE EQUATIONS FOR THE ORDER PARAMETER FOR A 
FERROELECTRIC WITH THE LIFSHITZ INVARIANTS 

The Landau free energy functional reads 

$ =Jdz {-r{rf, + r^l) + 711(77? + r^l)^ + U2{rilril)] + 
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where rj^ and are the order parameter components; we consider only one-dimensional configurations; parameters 
r, fj, ui, U2 and 7 are the Landau free energy expansion coefficient. Introducing the amplitude and phase variables 

7?! = pcosip,T]2 = psmip, 
we obtain the following expression for the Landau free energy: 

$ = / dz{— rp^ + up"^ + wp^ {1 + COB nip) 
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Varying the free energy we obtain the equilibrium equations Q 



rp + 2up^ + '^wp^ 1(1 + C0S7i(y3) +7p(^)2 - 
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nS^tp ^ dp dip dp n „ . 

+ ^^P^~^ '^PTT + o^^ smmp = 0. 

Here n is an integer number describing the system symmetry. Now we introduce the dimensionless variables: 

_ rr PF dR _dR FF 
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Then the equations ^ and (0]) take the form 

R" - R^ + (1 - </3'2 + T^')-R - 
R''-^K{cosnp+l) = 0, 



+ 2— - — T + i?"-2^sinn(p = 0. 
i? R 
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III. APPROXIMATE ANALYTIC SOLUTION 



Let us begin from the limit K = Q. Then the equation ^ can be solved: 

^ =^^=i?^+2' 



(8) 



where Co = (0) — ^] ^ (0)^ is the integration constant, which is determined by the initial conditions "0 (0) and 
R (0) . Now we can substitute this expression for p' into the equation ^ . We obtain in result a closed equation for 
the amplitude function R : 



i?"-i?3 + i?(l + ^)-^=0. 
This equation can be interpreted as a dynamics equation with the effective potential: 

(1 + -)- — + 



(9) 
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An interesting feature of this potential is its dependence on the initial conditions via the constant Cg . There exists 
a domain of parameters, where the potential has a minimum and, therefore, an oscillating solution for the amplitude 
function can take place. A condition of the maximum and minimum points merging into the inflection point with the 
horizontal derivative U' {R) — U" (R) — gives the equation 

y6 ^ ^2T-^ + 48T2 + 64 - 432C^ = 

However, a presence of a minimum is necessary but not sufficient condition: an oscillating solution will not exist if 
the initial point is situated outside the potential well. 

A border Hne for the domain, where oscillating solutions can exist is depicted here: 



FIG.l. A border line for the domain, where oscillating solutions can exist. 
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If we have a solution for the ampHtude function R{£,), the phase function (p{^) is given by the equation ([8]). The 
phase function consists from two terms: a slow contribution stemming from the second term in ^ and more or less 
diffused jumps due to the first term in this equa tion. It is important that the jump value equals exactly to tt : 

The first integral of [6] reads || = y/2{E~U), so = -j0=. 

Notice that the jump value / calculated in the vicinity of the R{£,) minimum equals to 

= 2f^'= — P — = 2arccos|#i| TT, (11) 

where Re 3> i?o- We have used the first integral of ^ above. Only leading terms of U{R) were taken into account. 
Note that the model ([2]), which can be reduced to the sine-Gordon model in the limit of i? = const ^ admits solutions 
with jumps exactly equal to ±7r (topological charge However, we have obtained a similar result within the 
approximation K = 0, i.e. neglecting the cosine term in ^\ As it is seen from the formulae flTl these phase jumps 
appear due to the excursion of the amplitude function near the singularity ^ . 

The equation ([9]) can be solved analytically, but we will present below numerical solutions both for K — and 



IV. NUMERICAL SOLUTION FOR K = 



Some results of calculation of the amplitude and phase functions spatial dependence for the case of if = are 
presented below. Figures 2-5 and 6-9 differ only by initial values of the amplitude function: in the case of figures 2 
- 5 we have small initial amplitude function and, therefore, oscillations of the amplitude function are near to harmonic 
ones, while in the case of figures 6-9 the initial amplitude function is near to the apex of the hump and, therefore, 
we have a train of solitons. 



FIG.2. Effective potential forn = 4, if = 0, T = 1, i?(0) = 0.3, R'{0) = 0, <^(0) = 0, <f'{0) = 0.3. 

Vertical dash marks R{0) value. 




FIG. 3. Spatial dependence of the amplitude function for n = A, K = 0, T = 1, i?(0) = 0.3, R'{0) = 0, <^(0) = 0, 

ip'{0) = 0.3 



FIG. 4. Spatial dependence of the phase function for n = 4, if = 0, T = 1, ii(0) = 0.3, i?'(0) = 0, ^'(0) = 

(a) thin Hne <f'{0) = 0.3 
(b) heavy Hne ip'{0) = 0.7. 




FIG. 5. Amplitude-phase polar diagram for n = 4, isT = 0, T = 1, i?(0) = 0.3, R'{0) = 0, ip{0) = 

(a) thin line <^'(0) = 0.3 
(b) heavy hne (f'{0) = 0.7. 
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FIG. 6. Effective potential for n = 4, K = 0, T = 1, R{0) = 1.099, i?'(0) = 0, <f{0) = 0, ip'{0) = 0.3 

Vertical dash marks -R(O) value. 
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FIG. 7. Spatial dependence of the amplitude function for n = 4, ii" = 0, T = 1, R{0) = 1.099, R'{0) = 0, ^(0) = 0, 

ip'{0) = 0.3 




FIG. 8. Spatial dependence of the phase function for n = 4, ii' = 0, T = 1, i?(0) = 1.099, R'{0) = 0, ip{0) = 

(a) thin hne <f'{0) = 0.3 
(b) heavy hne ip'{0) = 0.7. 
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FIG. 9. Amplitude-phase polar diagram for n = 4, ii" = 0, T = 1, i?(0) = 1.099, R'{Qi) = 0, <p(0) = 

(a) thin line !^'(0) = 0.3 
(b) heavy hne <p'(0) = 0.7. 
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V. NUMERICAL SOLUTION FOR K 

Results of the numerical solution of our equations in the case of K ^ are presented in the figures 10 - 15. The 
following important distinctions from the case of K = must be mentioned: 

(i) The periodicity of the spatial dependence is broken. 

(ii) The frequency and amplitude modulation can be seen in Figs. 10 and 11. 

(iii) A direction of the phase function staircase spatial dependence may change to the opposite after some number 
of steps. The number of steps between neighboring direction changes is random (see FIG. 15). 

(iv) There exists a parameter range, where the trajectories in the polar diagram show closed periodic movement: 
a synchronization with the periodic contribution of the sine of the monotonously increasing component of the phase 
function takes place. 




20 40 60 80 100 



FIG. 10. Spatial dependence of the amplitude function for n = 4, K = 1.6, T = 1, R{0) = 0.3, R'{0) = 0, ip{0) = 0, 

(^'(0) = 0.3 




FIG. 11. Spatial dependence of the amplitude function for n = A, K = 1.6, T = 1, R{0) = 0.3, R'{0) = 0, <p(0) = 0, 

<f'{0) = 0.75: thin line 
Spatial dependence of the parameter (see below) in arbitrary units: heavy hne. 




FIG. 12. Spatial dependence of the phase function for n = 4, K = 1.6, T = 1, i?(0) = 0.3, i?'(0) = 0, ip{0) = 



(a) thin line (f'{0) = 0.3 
(b) heavy hne <^'(0) = 0.75. 
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FIG. 13. AmpHtude-phase polar diagram for n = 4, K = 1.6, T = 1, R{0) = 0.3, R'{0) = 0, <p(0) = 0, <f'{0) = 0.3 




FIG. 14. Amplitude-phase polar diagram for n = 4, if = 1.6, T = 1, R{0) = 0.3, R'{0) = 0, ip{0) = 0,(p'{0) = 0.75 
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VI. DISCUSSION 

Thus, the main distinctions of the solution in the case of the non- vanishing anisotropy parameter K are a random 
change of the jumps direction (see FIG. 12) and an ampHtude modulation of the amplitude function (see FIG. 11). 
Just these random direction changes lead to the tangled trajectory in polar coordinates presented in FIG 14. 

In the case of if = we introduced the integration constant Co ([8]). Let us introduce the function: 



C(0 = 



Direct differentiation shows that this function reduces to a constant in the case of vanishing anisotropy parameter 



K: 



dC dR 



^ 2 



R^{-R"-^Ksmrnp) = -R^Ksiurvp. 



(12) 




FIG. 15. Spatial dependence of the parameter C for n = 4, K ^ 1.6, T = 1, R{Q) = 0.3, R'{0) = 0, 

ip{0) = 0,ip'{0) = 0.75 
(a) thin line C 
(b) heavy line C. 

Comparison of FIG. 15 and FIG. 12 confirms that changes of the direction take place in the points where C(^) 
function changes its sign. 



VII. CONCLUSION 



We have shown in this paper that the constant amplitude approximation gives a poor description of the real picture 
of the spatial evolution of the amplitude and phase functions for the model of the incommensurate ferroelectric with 
the Lifshitz invariant. 
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